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Abstract 

We generalize the analytical expressions for the two-loop leading-log neutral 
Higgs boson masses and mixing angles to the case of general left- and right- 
handed soft supersymmetry breaking stop and sbottom masses and left- 
right mixing mass parameters (tuq, mu,mD, A t , Ab). This generalization is 
essential for the computation of Higgs masses and couplings in the presence 
of light stops. At high scales we use the minimal supersymmetric standard 
model effective potential, while at low scales we consider the two-Higgs 
doublet model (renormalization group improved) effective potential, with 
general matching conditions at the thresholds where the squarks decouple. 
We define physical (pole) masses for the top-quark, by including QCD self- 
energies, and for the neutral Higgs bosons, by including the leading one-loop 
electroweak self-energies where the top/stop and bottom/sbottom sectors 
propagate. For tuq = rrtu = m^, and moderate left-right mixing mass 
parameters, for which the mass expansion in terms of renormalizable Higgs 
quartic couplings is reliable, we find excellent agreement with previously 
obtained results. 
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1 Introduction 



The use of effective potential methods has proved to be an elegant and simple way 
of incorporating all dominant top mass dependent one-loop radiative corrections to 
the Higgs masses and mixing angles in the minimal supersymmetric Standard Model 
(MSSM) |TJ, [|. Indeed, the one- loop effective potential computation reproduces with 
a good level of accuracy the complete one-loop radiative corrections to the Higgs 
masses |§ and mixing angles. However, relevant corrections are missing within this 
approximation. In particular, for the values of the top-quark mass preferred by the 
most recent experimental data ||, ||, the next-to-leading order Higgs mass values can 
differ from the leading order ones by 5 to 15 GeV. 

The two-loop corrections to the Higgs masses in the limit of large values of the 
CP-odd Higgs mass and stop mass parameters have been computed in different ap- 
proximations • In particular, a complete computation of the next-to-leading 
order effects on the lightest Higgs mass was presented in Ref. ||. It was subsequently 
realized that, when the leading- log renormalization group improved Higgs mass expres- 
sions |T(|-|TB|] are evaluated at the pole top-quark mass scale M t , they reproduce the 
next-to- leading order values with a high level of accuracy, for any value of tan j3 and 



of the stop mixing mass parameter |14 |. This means that, for this particular value of 



the renormalization scale, the genuine two-loop corrections are small, an observation 
analogous to the one already made in Ref. ||. The same holds if the renormalization 
scale is fixed at the on-shell top-quark mass, m t = m t (m t ), where m t is the running 
top-quark massQ. Based on the above observation about the choice of the renormal- 
ization scale, we presented in Ref. |TJ| analytical expressions for the two-loop leading 
order Higgs masses and couplings, which are valid for common values of the stop and 
sbottom supersymmetry breaking masses and moderate values of the squark left-right 
mixing mass parameters |T3|. Indeed, all next-to- leading order computations of the 
Higgs masses and couplings assumed the left-handed and right-handed stop mass pa- 
rameters to be equal. Moreover, they relied on an expansion of the effective potential 
which becomes reliable only for moderate values of the stop mixing parameter (see the 
discussion in Ref. [Q). It is useful to find a good analytical approximation, which 
works independently of the nature of the stop spectrum. The main objective of this 
work is to provide such an analytical approximation, based on an analysis of the domi- 
nant leading order effects in the effective potential computation of the neutral CP-even 
Higgs masses. 

We now present, for completeness, the neutral Higgs one-loop effective potential in 
the MSSM and define our notation and conventions. In a mass-independent renormal- 



ization scheme such as MS or DR [15|, the one-loop effective potential of the MSSM, 



as a function of the neutral components H 1 and H 2 of the two-Higgs doublets, is given 



lr The relation between the running mass mt and pole mass M t , for the top-quark, taking into 
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account one-loop QCD corrections, is: M t = mt(M t ) 1 ' 
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m~ and m 2 ~ are the mass parameters of the left-handed and right- 
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handed stops, while denotes the left-right stop mixing mass parameter. These 
mass parameters are explicitly given by 
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where the terms depending on the gauge couplings come from the effective potential D- 
terms. Similarly, in Eq. (|1.4j) m~ and m~ are the mass parameters of left-handed and 
right-handed sbottoms, and rrix b denotes the left -right sbottom mixing mass. They 
are given by 
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The mass parameters rriQ, rriu, mo, A t and At,, in (|1.5| ) and ( |TT6|) , are the soft breaking 
masses for the left-handed stop/sbottom doublets, right-handed stop and sbottom sin- 
glets, top and bottom trilinear couplings, respectively, while fi is the supersymmetric 
Higgs mass parameter. 

The stop and sbottom mass eigenvalues are given by 
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The content of this paper is as follows: In section 2 we re-analyse the effective 
potential computation of the neutral Higgs masses, putting special emphasis in the stop 
decoupling effects. For simplicity of presentation, we shall concentrate there on the case 
of large values of the ratio of the two Higgs vacuum expectation values, tan f3, and equal 
soft supersymmetry breaking parameters for the left- and right-handed stops, ignoring 
momentarily the bottom-sbottom contributions. The two-loop leading-log expression 
for the Higgs masses will be first presented within this framework. In section 3 we shall 
generalize the results of section 2 to the case of non-degenerate squark mass parameters 
and arbitrary values of tan f3 and the squark mixing mass parameters. The inclusion 
of the bottom- and sbottom-dependent corrections will be discussed. We also define, 
from the running masses of neutral Higgs bosons obtained from the effective potential, 
the corresponding pole masses by the inclusion of self-energies where we keep the 
leading contributions coming from the top/stop and bottom/sbottom propagation. We 
reserve section 4 for our conclusions, and appendices A and B for the explicit analytical 
expressions of the different self-energies and the gaugino/Higgsino corrections to the 
Higgs mass matrix elements, respectively. 



2 The limiting case of large tan j3 and rriQ = mu 

In this section we study in detail the simple limiting case of the MSSM with very large 
values of tan (3. This limiting value is more than an academic exercise since it provides, 
for a given set of the other free parameters of the theory, the absolute upper bound 
on the mass of the lightest CP-even Higgs boson and is therefore useful for practical 
applications. We shall set here the //-parameter to zero and neglect the bottom Yukawa 
coupling^ hb- Furthermore, we shall assume the case tuq = mu = M to simplify the 
analysis. All these constraints will be relaxed in the next section, where the general 
case will be studied. 

A good description of the effective potential is given by just taking Hi = and 
H 2 = H. For scales Q > M it can be written as 

V = K + m 2 \H\ 2 + -\\H\ 4 + V {1) (2.1) 
2 

where A cx M 4 is the one-loop vacuum energy induced by supersymmetry breaking, A 
is fixed by supersymmetry to 

x = \{g 2 + g' 2 ), (2.2) 

and is the one-loop effective potential (|1 . 1|) where only stop and top contribu- 
tions are taken into account. In the limiting case we are considering, and neglecting 

2 This is an excellent approximation even for moderately large values of tan/3. 



O (g 4 , g 2 g' 2 , g' A ) contributions in radiative corrections, the masses contributing to 
can be written as: 
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and the /3-functions of the tilded parameters in 
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For scales Q < M we have to properly take into account the process of stop decou- 
pling. In a mass-independent renormalization scheme [15] this can be explicitly done 
as follows^. We expand the mj 12 contribution to ( |2.1| ) around log(M 2 /Q 2 ) as: 
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and absorb the renormalizable terms in ( |2.5| ) in a redefinition of the low energy param- 
eters: 
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where all expressions are expanded up to one loop. 

Using now ( |2.4j ) and (215) we obtain the one- loop /5-functions for the untilded 
parameters 
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Notice that ( |2.7| ) are the one-loop /5-functions of the Standard Model (as they should) 
to the order of approximation we are working. In the same way H = H and h t = h t , 
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since neither the anomalous dimension nor the top-quark Yukawa coupling are affected, 
in our approximation, when passing through the threshold scale M. 



as 



Using now ( |2.1| ), ( |2.5| ) and Q2.6| ) we can cast the effective potential for scales Q < M 
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where G(\H\) = 0(\H\ ), which contains the higher dimensional terms arising from 
( p.5|) , is scale independent and can therefore be 'frozen' at the scale Q = M. It is given 
by 
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The effective potential ( |2.8j ) is a function of a; 
ij} = ImH. Imposing the condition of minimum at (0, i/j) 
v(m t ) = 174.1 GeV 



H\ 2 = 



4> 2 + ip 2 , where = Reif , 
v , 0), where we are taking 
14| , one obtains a massless scalar (neutral Goldstone boson) 
along the direction ImH. The condition of minimum amounts to dV(x)/dx = and 
can be written, at Q 2 
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where we are neglecting two-loop corrections. 

Imposing now condition ( |2.10|) one gets for the Higgs mass 
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The usual procedure is to neglect the term d 2 G/dx 2 in (|2.11|) , in which case the 
usual expression for the Higgs mass squared in the leading order approximation, 2\v 2 , 
is obtained. We shall now take into account the contribution from the G-function and 
see how it modifies the final result. We shall first analyse the case of zero mixing. 



2.1 The case of zero mixing 

In this case A t = 0, and the expression for G in (|2.9|) simplifies a lot. Using ( |2.6|) and 
one can, to one loop, expand A as 
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where all couplings are considered at the scale M (remember that the G-term was 
frozen at that scale) and we have neglected 0(g 4 ) terms. 

Putting (|2.13|) and (|2.14|) together, the effective coupling (|2.12|) can be written as 
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where all couplings on the right-hand side are considered at the scale M. 

Comparison of ( |2.13| ) with ( |2.15| ) shows that the effect of including the G-term in 
( |2.8|) is, apart from the last (mixed Yukawa-gauge) term in ( 2.1 5| ) which is numerically 
unimportant, to resum the series 
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from the expansion parameter log(M 2 /m 2 ) to the physical one log((M 2 + m 2 )/m 2 ), 
which goes to zero in the supersymmetric limit M — > 0, and makes the radiative 
corrections vanish in that limit. 

In other words the resummation from the G-function allows us to replace the ex- 
pansion ( |2.13j ) by the expansion 



A eff (m t ) = X cS {M s ) - P\ eff (M s )t 
where the expansion parameter is 



and 



t 



Ml 



, M 2 



mi 



M 2 + m 2 . 



(2.16) 

(2.17) 
(2.18) 



2.2 The case A t ^ 
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In the absence of the G-contribution in ( |2.12| ), Eq. ( |2.19| ) gives the usual one-loop 
leading-log contribution to the Higgs mass. Introducing now the G-function, one can 
write ( [2.12j ) as: 
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where the 'threshold' contribution to the coupling is given by 
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where the ellipsis denotes the contribution from higher- dimensional terms in m 2 /M| 
and A 2 /Ml 

We can see that the effect of taking into account the G-function in (|2.8| ) is twofold: 



It makes a resummation from log(M 2 /m 2 ) to the physical expansion parameter 
log(M|/m 2 ) in the one-loop radiative correction, as well as resummation from 
A 2 t /M 2 to A 2 /Ml in the threshold term of ( g!9|) . 



It generalizes the threshold correction of Eq. ( |2.19| ) to include higher order effects 
in powers of A 2 /M 2 S in Eq. ( |2T2TD . 



Therefore the net effect of the G-function is to change the expansion (gl9|) into 

A eff (m t ) = A cff (M 5 ) - (3 XcB (M s )t (2.23) 

where 
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2.3 Two-loop leading-log expansion 

As we have seen in Ref. [[L4j one can obtain the two-loop leading-log correction by 
expanding the parameter A to order (log(M 2 /m 2 )) 2 , as 



M 2 1 , / M 2N 
A(m t ) = A(M) - /3 A (M) log ^ + -/3 A log 
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where, since we are now considering two-loop corrections, we have evaluated the quartic 
coupling at the on-shell top quark mass scale and the prime denotes derivative with 
respect to log(Q 2 )- In the previous subsection we have seen that the effect of including 
the G-function in the effective potential is to resum M 2 to M| in the first two terms 



of ( 2.25 ) and replace the threshold correction in the first one by A t hA e fr in ( 2.21|) (in 



Ref . |L4| , this resummation effect was assumed to be true to assure a proper behaviour 
of the radiative corrections for M = 0). To prove resummation in the last term of 
( p.25|) one would need to use the whole two-loop effective potential in the MSSM. In 
the absence of such a calculation we shall assume that this happens, since we have 
already proved that t in ( p. 17 ) is the physical expansion parameter in the one-loop 
calculation, and in addition the numerical relevance of the resummation in the two- 
loop corrections is expected to be tiny. We shall therefore consider as a starting point 
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where A e ff(M s ) is given in ( |2.24| ) and t and M s are obtained from Eqs. ( |2.17| ) and 
( B - Defining (3 XcS = a x 

eff A e g- + b\ aS , it follows that 



Ks(m t ) = Ks(M s ) (1 - a XcH (m t ) t) - b XcB (m t ) t (1 - a XcS (m t ) t) - ^(3' XaH (m t )t 2 . 

(2.27) 

From the renormalization group equations of the quartic coupling, it follows that the 
term 1 — a Xcff (m t ) t = ^ _4 (m ; ), where £, the Higgs field anomalous dimension, is: 



Km) = 1 + ^ h t( m t) L 



Recalling that v 2 (M s ) = v 2 {m t ) £ 2 (m t ), from Eq. (|2.11|) we get 

m 2 h (m t ) = m 2 h (M s ) £~ 2 (m t ) + A rad m 2 h (m t ). 



(2.28) 



(2.29) 



In the above, 

m 2 (M s ) = 2\ cS (M s ) v 2 (M s ), (2.30) 

X e s(Ms) is given in Eq. ( |2.24| ) with all couplings and masses evaluated at the scale Ms, 
and A ra d m 2 (rn t ) is given by 
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where all couplings in ( |2.31| ) are evaluated at the scale Q 2 = m 2 , a 3 is the strong gauge 



coupling, and we have used (3h t in the evaluation of ( |2.31|) Observe that, if we 

replace A t hA e ff by its expansion in powers of At/ Ms, Eq. (|2.22|) , we neglect the small 
terms depending on the weak gauge couplings, and we re-express the values of the 
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couplings at M s by their expressions at m t using the appropriate 7- and /^-functions, 
we obtain 
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Equation (|2.32|) is equivalent to Eq. (9) of Ref. |14j] in the large tan (3 regime. 

Although we have used, as simplifying hypothesis, the case where tan/3 3> 1 and 
\i — 0, all the results are also valid for the case iua ~ Ms and any value of tan/3 and 
fi. The only change in the final results is that 



(g 2 + g' 2 ) 



(g 2 + g' 2 ) cos 2 2/3 
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in A(M) and A eff (M s ). 

Finally, we want to conclude this section with a comment about the physical in- 
terpretation of the decomposition fl2.29| ). The term m^Ms) comes from the scale 
independent part of the MSSM effective potential frozen at the scale Ms, where we 
have already subtracted the contribution evolving with log(M|/m 2 ). Since this term 
is scale independent, it is evolved to the scale m t with the corresponding power of the 
anomalous dimension of the Higgs field. On the contrary, the term A rac jm 2 , which 
arises from renormalizable terms after the stops are decoupled at the high scale, is 
computed at the low scale m t . Had we considered the whole MSSM effective potential 
at the scale m t 0, we would have been neglecting the stop decoupling at the scale Ms- 
This is shown in Fig. 1 (dotted lines) where we plot as a function of Ms for a pole 
top-quark mass M t = 175 GeV, vanishing mixing, A t = fi = 0, and large (tan/3 = 15) 
and small [infrared (IR) fixed point solution: sin/3 ~ (200 GeV/M t )] values of tan/3. 
Had we evolved the Higgs mass obtained from the whole MSSM effective potential 
(including the logarithmic terms) at M s , to the scale frit with the anomalous dimen- 
sion factor t;~ 2 (m t ), we would have made an error associated with the non-exact scale 
invariance of the effective potential in the one-loop approximation, as was observed 
This is shown in the dashed lines of Fig 
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1 and was also noticed in Ref 
In fact the latter procedure would lead to an expression of A rad ml, where the second 
term inside the square brackets has an extra factor of 2 [see Eq. ( |2.32j) 1. The solid lines 
in Fig. 1 show the corresponding value for the Higgs mass, while considering Eq. ( |2.29| ) 
with the whole expression for A t hAA cff from Eq. ( p. 21 



The dependence of ( |2.29p on the mixing is shown in Fig. 2 (solid lines) where we 
plot rrih as a function of A t , for /i = 0, the same values of tan/3 as in Fig. 1 and 
Ms = 1 TeV. For comparison we also plot (dashed lines) the corresponding mass, 
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using the approximate expression for the threshold contribution to the Higgs quartic 
couplings, Eq. ([2.32 ). We see that this approximation, which was used in [14]j, is very 
good up to the maximum of the curve, which means that the absolute upper bound on 
the mass of the lightest Higgs boson previously obtained remains unchanged for any 
value of the mixing. Of course for very large values of the mixing there is a departure 
between the two curves. We have deliberately omitted from Figs. 1 and 2 the small 
D-term contributions in Eq. ( |2.21| ) to compare with previous results which did not 
make use of them 0, [14]]. They will be included in the general analysis of section 3. 



3 The general case 

Having understood the simplified case explained in section 2, we can proceed with 
the general case. Let us first focus on the behaviour of the renormalized Higgs quar- 
tic couplings for low values of the CP-odd Higgs mass and values of the mass pa- 
rameter vhq different from mu- In this case, the effective potential at scales above 
Max(mQ, mu, m^) is obtained through the general expression, Eq. ( |1.1| ). Expanding 
the effective potential in powers of the Higgs fields, it is easy to identify the form of 
the effective quartic couplings at the scale mq, Aj(mg), in a way completely analogous 
to what we have done for the case mq = mu in section 2. Assuming for definiteness 
that mq > mu > mo, 

X t {m Q ) = \irnu) + (3? U (mu) t QU + (j3? u )' ^f- (3.1) 

[see e.g. the second equality in Eq. ( |2.26| )1 where tQu = log(mg/m^), fa denotes the 
/3-function of the Higgs quartic couplings, fa = d\i/d\og(Q 2 ) and the superscript QU 
denotes its expression at energy scales between mq and mu- In general, we can write, 

fa = diXi + b { , fa ~ Oibi + 6J, (3.2) 

where we have only kept the dominant, Yukawa-coupling-dependent contributions to 
fa. We omit the scale dependence of ft[ since it is a higher-order effect. Hence, 

K.{mu) = K.{m Q ) (l - di t QU ) - bf u {mu) t QU (l - Oj t QU ) - (f3® U ) (3.3) 

The coefficients Oj are linear combinations of the anomalous dimensions of the Higgs 
fields Hi, which are independent of the squark fields. The same procedure as used above 
may be used to connect the value of the quartic couplings at mo with their values at 
mu, and finally the values at mu with their values at m t . Keeping only the dominant 
terms, we obtain 



QU\' tQU 
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( ' U (m) t QU 1 - ai(m t ) (t Q + ti 



(t Q + t v 



- b. 
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- b. 



> tr 



1 D (m t ) t UD [l - a t (fn t ) (t v + t D )} - (/3f D ) (t v + 



~ t 

H m t) (l - ai(m t ) t D ) - & -y-, 



(3.4) 



where /3j without any superscript denotes the /^-function values once the two stops and 
the two sbottoms are decoupled, the superscript XY, with X, Y — Q, U, D, denotes 
the functional form of the /^-functions at scales between mx and my, txy — tx — ty 
and tx = \og(m 2 x /rn 2 ) . Similar expressions are obtained for a different hierarchy of the 
squark mass parameters, with the only difference that, in the case mu > mq and/or 
mo > rriQ, stops and sbottoms should be decoupled at different scales. For simplicity 
of presentation we shall first discuss the result for the Higgs mass matrix elements in 
the case under study and we shall present below the result in the most general case. 

The contribution of the quartic couplings to the Higgs mass matrix elements is then 
given by 



M 2 2 
M\ x 
Mf 2 



+ A 6 cos 2 (3 + X 7 sin 2 0\ 
- 2A 6 cos /3 sin (3 + A 5 sin 2 (3\ + ™ 2 
2v 2 [A 2 sin 2 (3 + 2 A 7 cos (3 sin (3 + A 5 cos ; 



2t> 2 [sin (3 cos/5(A 3 + A 4 
2w 2 [A!Cos 2 /?-' 



„2 



i 2 A sin (3 cos (3 
m 2 A sin 2 (3 
2 f3] + m\ cos 2 (3, 



(3.5) 



where all terms should be evaluated at the same scale and we have also included the 
dependence on the CP-odd Higgs mass [O. In the above, v i = v cos (3 and v 2 = v sin (3 



are the Hi and H 2 vacuum expectation values, respectively. 

Equation ( |3.4|) has a clear interpretation: The factor X^rnq) in the first term con- 
tains the tree level terms and all finite contributions to the quartic couplings, arising 
from the existence of squark mixing and the fact that mn ^ vtlq ^ mu. The factor 
involving <2j in the first term contains exactly the terms necessary to rescale the Higgs 
mass matrix elements by the appropriate anomalous dimension factors from the scale 
mq to the scale m t , together with the ones necessary to re-express the vacuum expec- 
tation values of the Higgs fields appearing in Eq. ( |3.5|) at the scale frit, in terms of their 
values at the scale rriQ (for the complete expression of the /^-functions of the Higgs 
quartic couplings, see for example Ref. [0)- For instance, singling out this contribu- 
tion of the A2 coupling to the matrix element M| 2 , which we shall denote by K22, we 
obtain 







1 - a 2 log 






\m 2 ) _ 



2 A 2 (m Q ) 
2 A 2 (m Q ) i 2 A im t ) v\(m t ) 



v 2 (m t ) 



(3.6) 



= 2 A 2 (m Q ) vl(m Q ) £ 2 2 {m t ), 

where £ 2 is the anomalous dimension factor of the Higgs fields H 2 . The contribution 
of the other couplings to the Higgs mass matrix elements present similar properties. 
Thus, both the tree level term and all finite terms leading to the non-trivial matching of 
the quartic couplings at the scale mg may be treated in the same way as the terms pro- 
ceeding from the higher- dimensional operator contributions to the Higgs mass matrix 
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elements. That is, they may be frozen at the scale tuq and rescaled with the appro- 
priate anomalous dimension factors to obtain their expressions at low energies. This 
generalizes the result obtained in section 2 for the case of degenerate squark masses, 
itlq = mu- 

The following terms in Eq. (|3.4|) are the ones which would be obtained even in 
the presence of trivial matching conditions for the quartic Higgs couplings and, as has 
been clearly explained in the previous section, are associated with the scale-dependent 
contributions to the effective potential. The dominant leading-log contribution to the 
quartic couplings proceeds from the terms in the /3-function proportional to the fourth 
power of the top and bottom quark Yukawa couplings, which are given by 




A,(m t ) - K(m Q ) = log + log , (3.7) 



where b\ = —3hf/8ir 2 and b\ = —3h^/8TT 2 , and h t are the bottom and top quark 
Yukawa couplings, and X = D, U for i = 1, 2, respectively. Although 6 3 and 64 also 
present a quartic dependence on the Yukawa couplings, such dependence is absent from 
63 + 64. As may be easily proved using Eq. ( |3.5|) , and following the same procedure as 
in section 2, in the case of no mixing the contribution of the higher-order operator to 
the '22' ('11') Higgs mass matrix elements allows us to replace the factors wig and m\j 
(tUq and m 2 D ) in the leading order expressions by mQ + m 2 and mfj + m 2 (rag + ml and 
m 2 D + m|). The same occurs with the D-term contributions proportional to the square 
of products of the weak couplings and the Yukawa couplings. It is hence convenient to 
define 



t Q = log (^^=2 — J and tu = lo S 

fS = los (^\ and fD =lo g (^] (3.8) 



while tQu = to — tu and tqn = tR — tp. Observe that in the denominators of the 
expressions for t D and tR we have written instead of m t since we know that, upon 
consideration of the whole supersymmetric contributions to the physical masses, ra- 
diative corrections should vanish in the supersymmetric limit; we should thus use 
expansion parameters tq and tp, which vanish in the limit tuq, vtld — > 0. This change 
has negligible effects on the Higgs mass computation. 

In the general case, the way to proceed to obtain the Higgs mass matrix elements 
at the scale m t is the following: the CP-even Higgs mass matrix elements may be 
decomposed in three terms, namely: 

Ml(m t ) = Ml{m t ) + (^.(M st )) T (eKm t )ej(m t ))" 1 

+ (^(M s6 )) T (e[(m t )g(m t )) _1 , (3.9) 

where £i(m t ) denote the anomalous dimension factors, 
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where for convenience we define, 



1 + 



3hl 
32tt 2 



ft 
' i 



t b 1 = Max(tQ,t D ) t\ 



Max(t%,tu), 



and also, for later use, 



Min(tg,t D ) 



Min(tg,% 



(3.10) 



(3.11) 



(3.12) 



In the above, ^Mf 3 -(M st )J~ ^A^f^Mj^j-J is the contribution to the mass matrix ele- 



ments coming from the terms frozen at the scale M st (M^), where the stops (sbottoms) 
are decoupled, with M 2 t = Maxfrn^+mj, m^+m 2 ) (M 2 b = Max(mg + m 2 , m 2 D + m^)), 

and are obtained from the mass matrix elements M 2 ^, Eq. ( |3.5|) , by considering the 
one- and two-loop leading logarithm contributions in the renormalizable Higgs quartic 
coupling expressions. For example, in the case mq > mjj > mjj, these contributions 
to the quartic couplings Aj are given by 



Xi = A; - (\i{m Q ) - Af 0C (m Q )) (l - ai(m t ) t Q ) 



(3.13) 



where the A, are given in Eq. A*/ 66 = A*/ 66 = - (Af 56 + Af 36 ) = (g 2 + g' 2 ) /4, and, 

in order to simplify the presentation, we include all D-term contributions in the defi- 
nition of the quartic couplings \. Observe that, since M^iyrit) contains only the one- 
and two-loop leading logarithm expressions independent of the mixing mass terms, the 
quartic couplings A5, A6 and A7 give no contribution to M^imt). Replacing the quartic 
coupling /3-functions by their dominant Yukawa coupling dependence, we obtain, in 
the general case 



Ai 



+ 



+ 



g 2 + g' 2 



2 h t 



8tt 



-h 2 t b 

• a b l i 



16n 
8tt 2 hb h 



r - 
1 + 



16tt 2 



-h 2 b 
2 b 



^h 2 ~8g 2 3 



t\ + 4 



16vr 2 



-h 2 b 
2 b 



^h 2 -8gj)t 



D 



(3.14) 



+ 



+ 




8 g\ (t\ + 4 



(3.15) 



A3 + A4 



g 2 + g' 2 



D 



-h 2 t b 



16n' 



16n' 



-h 2 f 



(3.16) 
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where we have already performed the top and bottom mass resummations, leading to 
the change tx — > tx, with X = Q,U, D. In the above, h t and hb denote the top and 
bottom quark Yukawa couplings at the scale m t and #3 is the strong gauge coupling 
at the same scale. The Af terms are the additional leading-log D-term contributions 
appearing through the /3-functions of the quartic couplings, which contain additional 
terms proportional to the square of the product of weak gauge couplings and Yukawa 
couplings. Their two-loop leading-log contributions are very small and can be ignored. 
Interestingly enough, once these terms are considered together with the terms written 
explicitly in Eqs. (ETT^-flOj), one recovers the expressions for the D-terms first found 
in Ref. ||17|1 . Defining = ^A^ + (A^)^, for vtlq > mc, we obtain 



1 



16tt 



= ~^' 2 hlt QD (3.17) 



while for mo > tjiq we get 



K) 



D 



A D - 




K t QD (3.18) 



Analogously, for mq > mu we obtain 



8tt z 



QU 



*%)u = -^9' 2 hlt QU (3.19) 



while for mu > mQ we get 



A D 




A 3 L = -^ U 2 + V ^ tQU (3.20) 



32tt 2 



Finally notice that tan/? is fixed at the scale tua, for rriA < mt, while for > m t , 
tan p is given by 



tan/3(m t ) = tan/?(m^) 



(3.21) 



For the case in which the CP-odd Higgs mass tua is lower than Ms = Max (M st , M s b), 
we should have decoupled (strictly speaking) the heavy Higgs doublet at the scale m^, 
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and denned an effective quartic coupling for the light Higgs as A(m^) = m h (mA)/2v 2 , 
the low energy value of it being obtained by the running of the Standard Model renor- 
malization group equations from the scale rriA to m 4 [T3|. For simplicity we have 
ignored, in our analytical approximation, the effect of the heavy Higgs doublet decou- 
pling at the intermediate scale. We partially compensate this effect by relating the 
value of tan (3 at the scale m t with its corresponding value at the scale uia through its 
renormalization-group running, Eq. ( |3.21| ). 

The Higgs mass matrix elements at the scales M st , M s b may be inferred from the 
second derivative of the one-loop effective potential, Eq. at the minimum values 

for the Higgs fields. Consequently, they can be obtained from the expressions given 
in Refs. 0, [17|, where all parameters should be assumed to be given at the scale at 
which the matrix element is evaluated, i.e. M st or M sb . As we have shown above, the 
dominant D-term contributions to the Higgs masses are already taken into account 
in the one-loop leading logarithmic expressions, much as it happens in the case of 
irtq = mjj, where the dominant D-term contribution comes from the first two terms 
in Eq. fl2.2U| ). The D-term contribution of A t hA c g becomes relevant only for very large 
mixing and gives a very small correction to the top quark Yukawa coupling effect. For 
completeness, we shall include all 0(g 2 h 2 , g' 2 h 2 ) D-term contributions, where q = t,b. 

In order to obtain the expressions for (yVf^)~ , i^M 2 ^- one must subtract from the 
expressions of the matrix elements at the scale M st , M sb , respectively, the contributions 
coming from the term dependent on the CP-odd Higgs mass and from the leading order 
logarithmic expressions. This procedure leaves in the matrix elements all the terms 
that should be frozen at the scales M st , M s &, including all terms depending on the 
squark mixing mass parameters. Using the expressions given in Refs. [0, [T7|] , and 
Eqs. (gg)-(gTg), we obtain, 



K(M st )Y 



]ir 2 v 2 



A* U + (A' U )_ Af 2 + (A' 12 
A[ 2 +(A' 12 ) t Af 2 +(A 22 y 



(3.22) 



(M%{M sh ) 



]ir 2 v 2 



Al + (A' U )* A$ 2 + (A' 12 
A' 12 )' Af 2 + (A 22 )' 



(3.23) 



where 



+ 



ml 



cos 2 (3 



ml 



cos 2 (3 



log 



2 2 

m~ m~ 

bi h 



'm 2 Q + ml) {m 2 D + ml 



Af,(Ab — /itan/3) 



m- 



b 2 



, 2A b {A b - fj, tan P) . 
H o o log 



m- 



m- 



b 2 



(3.24) 



m. 



sin (3 



fi(—A t + [i cot (3) 



m- 



m- 



q(m 2 ~ , m 2 ~ x 
yv ti' ti> 



(3.25) 
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ra; 



A t (A t -//cot/3) 



sin 2 P 



ra- 



g(m 2 ~,m 2 ~), 



(3.26) 



*2 



ml 



//(— A b + //tan/3) 



'22 



COS 2 P 



m~ 



bi 



b 2 



g(m 



(3.27) 



m\ //(— A t + /i cot P) 



sin /3 



m- 



*2 



log 



ra- 



+ 



At (A t - //cot/3) 



ra- 



*2 



ra- 



m- 



m r 2 ) 



*2 



(3.28) 



//(— Ab + //tan/3) 
cos 2 /? m~ — m 2 



6i 



l>2 



//tan/3) 



f>2 



ra- 



ti 



ra- 



y(m 2 ~,m 2 ~) 



f>2 



(3.29) 



M 2 



Mf 



2m 2 J 1 b - ra 6 A b / 2 6 ] , (A' n )* = M 2 z m tl icot ptf 

-2m 2 J* + mtAtfj] , (A 22 ) 6 = -M 2 ra b //tan/3/ 2 ' 



m 2 cot — m t 



At cot P + // 



A 



2 6 A 6 tan/3 + // r 
-mjtan/J/! + m 6 / 2 



(3.30) 



?K, ™ 2 ) = 2 2 2 lQ S - 2 ' 

mf — m 2 m 2 



(3.31) 



A* 



ra- 



ti 



ra 



1 4 . 2 \ /ra 



ti 



*2- 



/I 2 o \ / rarrar 
+ sin 2 % log 2 fa % 

V2 3 / \mQ + mf / 

2 o / rarrar \ 

+ -sinolog' 1 " 1 



mjy + m 2 y 



(3.32) 
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+ 



m t 



A t — fi cot P 



ru- 



mp 



- sin 2 W 





fm 2 ~\ 


-\ log 







m Q - ™u 

i i 

m~ — m~ 
h t2 



g(m\,m\) 



(3.34) 



/a 



At — fi tan /? 



2 2 

m~ — m- 



b 2 





( m 2 ~ \ 


~ log 






\ m \) 







+ 



2 
3 



- — - sin 2 



m- 



m- 



-g{m. 



(3.35) 



In the above, all terms should be computed at the scale M st or M sb depending on 
whether they are associated to stop or sbottom contributions and all 0(g A , g 2 g' 2 , g' 4 ) 
terms are ignored. It is easy to show that, apart from very small terms of order 
M|/(m~ +m~ ), Eqs. ( |3.22| ), ( ^.23| ) vanish in the case of zero squark mixing. It is also 
straightforward to show that in the limit of large tan/3, run = rujj and mi, = u = 0, 

A{ 2 . md 3 (*< + (A,) 1 ) /16 .V reproves th e e X pl ion o f ^ gl ve„ ta 

Eq. ( |2.21| ). All parameter values in the expressions for the matrix elements .M 2 ^ can be 
expressed in terms of their values at the scale rh~ t by making use of the corresponding 
/3-and 7-functions. 

Having computed the renormalization group improved Higgs mass matrix elements 
at the scale rut, the neutral CP-even Higgs mass eigenvalues can be easily derived. 
They read 



m h(H) 



TtM 2 =f J(TtM 2 ) 2 - 4detX 2 



where 



TtM 2 = M 2 n + M 2 22 



detM 2 = M\ l M 2 22 



M 2 
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(3.36) 
(3.37) 



and the M.ij are the renormalized Higgs mass matrix elements at the scale m t , Eq. ( ^.9| ). 
From the matrix elements, the mixing angle a is also determined by [|13j : 



sin 2a 



2M 2 



12 



(TiM 2 ) - 4det A4 2 



(3.38) 
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cos 2a 



M\ 



Ml 



22 



(3.39) 



sj (TrA4 2 ) 2 



4detA4 2 



Concerning the running of the squark mass parameters, the dominant contribution 
comes from gluino-induced effects, which are absent in the case of heavy gluino par- 
ticles, as we are considering within this work. The remaining contributions are small 
and have a somewhat complicated dependence on the squark and Higgs spectrum. We 
shall ignore them within our approximation. We have further defined the light stop and 
sbottom masses as the values obtained using Eqs. (|1.7|) , (|1.8| ), while taking the running 
mass parameters at the scale m t and adding the QCD-dependent vacuum polarization 
effects. A more precise definition of the squark masses may be obtained by computing 
the squark effective potential and adding the full vacuum polarization contributions to 
the squark masses, much as we have done in the case of the Higgs bosons. We shall 
concentrate on this subject elsewhere. 

Figures 3 and 4 show the dependence of the Higgs mass m^ for varying values of 
itlq = A t = rriA and for a fixed value of the right-handed mass parameters my = 
mo — 1 TeV and my = my, = 100 GeV, respectively. The supersymmetric mass 
parameter fi has been set to zero. The solid lines represent the value of the Higgs 
mass by performing the renormalization group improvement of the effective potential 
method, as explained in this work. The dotted lines represent the values obtained 
from the one-loop effective potential, while ignoring the stop decoupling and taking all 
values to be given at the scale rn t 0. The dashed lines are the values obtained for 
the Higgs mass, while considering that the effective potential is scale-invariant, that 
is by taking the second derivatives at the scale Ms and rescaling them to the scale 
m t through the appropriate anomalous dimension factors (see e.g. Ref. Q)- F° r l° w 
values of my and mo, as those shown in Fig. 4, the last method becomes accurate for 
all values of mg < 600 GeV, while in the second method the departure from the proper 
renormalization group improved values is faster. Observe that the behaviour shown in 
Fig. 4 is very similar to the one shown in Fig. 1. Since the main purpose of Figs. 3 
and 4 is to compare the results in this paper with other different approaches, based 
on the one-loop MSSM effective potential, we have plotted in them the running Higgs 
mass at the scale m(, and turned the small D-terms threshold corrections off, i.e. we 
have put ■ = 0. In the following we shall turn these D-terms on and consider pole 
Higgs masses. 

The Higgs masses mh,u,A defined in Eq. Q3.36Q are all running masses obtained 
from the effective potential^, and evaluated at the top-quark mass scale. To compute 
the physical (propagator pole) masses M h ^ H> A one has to correct for the fact that the 
effective potential is defined at zero external momentum. In fact, the pole and running 
Higgs masses are related by (see e.g. Ref. [[|) 



5 In fact, as was noticed in Ref. [Q, the mass of the CP-odd Higgs, mj, turns out to be scale 
independent at one loop. 




(3.40) 



where ip = h, H, A and 




(3.41) 
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II (/ ,(g 2 ) being the renormalized self-energy of the corresponding Higgs boson. We have 
computed the Higgs self-energies, at the one- loop level, from the top/stop and bot- 
tom/sbottom sectors. The corresponding expressions can be found in Appendix A. 

Figures 5-7 show the variation of the pole Higgs mass Mh as a function of thq for 
fixed rriA and several fixed values of mjj and of the stop mixing parameter A t . Although 
in general, for a fixed moderate value of my and a fixed value of A t , the Higgs mass 
increases together with tuq, for large values of A t and moderate values of ray or for 
small values of my and moderate values of A t , situations can be observed for which the 
Higgs mass decreases with larger values of tuq. An understanding of this effect may be 



obtained by making use of the approximation of Ref. |T^| : although log(M S u SY /m t ) , 
with M| USY = (jn 2 ^ + m \^) /2, increases together with the squark-mass parameters, 
leading to larger values of the Higgs mass, the squark mixing contributions are maxi- 
mized for values of A^ ax ~ 2.4Msusy- Hence, for fixed values of A t a delicate balance 
between these two effects should be present in order to maximize the Higgs mass. 
Observe also that for small values of my, such as the ones analysed in Fig. 7, the 
configurations which maximize the lightest CP-even Higgs mass Mh correspond to sit- 
uations for which one of the stops may become light. Indeed, the curves are cut since 
we have introduced the (crude) experimental constraint mj 2 > 45 GeV on them. An 
important effect arising in the case of small values of one of the mass parameters, for 
instance my ~ M§, and tuq ^> mfj as shown in Fig. 7, is that for the same value 
of Msusy, the maximal value for the Higgs mass is always lower than in the case of 
rriQ = my. This is due to the fact that the requirement of having stop masses above 
the present experimental bound implies A t < nig ~ \/2Msusy- Hence, A t is always 
significantly lower than v4™ ax , which, as explained above, is the value that maximizes 
the Higgs mass for that particular value of the average stop mass scale Msusy- 

One could enquire about the stop and sbottom vacuum polarization contributions 
in the pole Higgs boson mass definitions. We have checked that these contributions 
do not give a significant effect in the determination of the neutral Higgs boson masses, 
unless one of the squarks becomes light (i.e. t and/or b) and its couplings to the Higgs 
fields (i.e. A t , At,, fi) are large. This behaviour is illustrated in Fig. 8 where we plot the 
neutral Higgs boson running (m h , m H , m^) masses (dotted lines) and pole (M h , M H , 
Ma) masses (solid lines) as functions of A t , for fixed values of the other supersymmetric 
parameters. We also plot the mass of the lightest stop mj 2 (dashed line). We see that 
the departure between the running masses and the pole masses occurs in all cases for 
large values of A t and small values of rnj 2 - 

The dependence of the Higgs mass on mo becomes relevant only for large values of 
tan/3. In Fig. 9 we show the variation of Mh as a function of mo for tan/3 = 60, fixed 
values of mQ, my and m^, and two values of ji: /i = 1 TeV (solid curves) and /i = 2 
TeV (dashed curves). The radiative corrections to M h induced by the bottom/sbottom 
propagation are always of negative sign and become only relevant for very large values 
of the /^-parameter. The behaviour with fi can be understood from the fact that the 
larger /i, the larger mixing in the sbottom sector and hence lighter sbottom masses 
are obtained. In this example, we have chosen mb(m t ) = 3 GeV, which corresponds to 
a bottom quark pole mass Mb ~ 5 GeV. All the curves are cut by the experimental 
constraint m-^ <; 45 GeV. Finally Fig. 9 also exhibits the dependence of the Higgs 
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mass on the parameter A t , its dependence with respect to A b being tiny. 



4 Conclusions 

We have presented a renormalization group improvement of the effective potential 
computation of the neutral Higgs masses in the MSSM. The method provides the first 
calculation of two-loop leading order corrections to the Higgs masses valid for any 
value of the soft supersymmetry breaking squark mass parameters, rriQ, mu, rn D , A t 
and Ab, the CP-odd mass ttia, the super symmetric Higgs mass \i and tan/?. This 
generalization is essential for the computation of the Higgs masses and mixing angles 
in the presence of light squarks. Our method uses explicit decoupling of stops and 
sbottoms at their corresponding mass scales, leaving threshold effects in the effective 
potential (and coupling constants) frozen at the decoupling scales and evolving, in the 
squared mass matrix, with the anomalous dimensions of the Higgs fields. 

The threshold effects achieve a complete matching of the effective potential for 
scales above and below the decoupling scales, and include all higher order (non- 
renormalizable) terms arising from the whole MSSM effective potential. The effect 
of considering non-renormalizable threshold effects in the effective potential is twofold: 
on the one hand it triggers resummations in the renormalization group expansion of the 
parameters, leading to 'physical' expansion parameters; on the other hand, it enables 
to consider the general case of arbitrary left-right squark mixing, as well as general 
left- and right-handed soft supersymmetry breaking squark masses. 

We have corrected the running neutral Higgs boson masses with one-loop self-energy 
diagrams, where top- and bottom-quarks, stops and sbottoms propagate, to define the 
corresponding pole masses. The numerical effect of polarizations is relevant only under 
special circumstances: light squarks and large mixing. 

We have analysed the general pattern of Higgs masses for general values of the 
supersymmetric parameters. We have found regions in the parameter space where the 
radiative corrections become large and negative. They are characterized by large values 
of the mixing-mass parameters, where the stability of the electroweak minimum can 
be endangered by the presence of charge and color breaking minima. Our results also 
allow the evaluation of the relevant radiatively corrected Higgs couplings through the 
corresponding value of the Higgs angle a [14[ . 



We have neglected, throughout the whole calculation, the possible contribution 
coming from light charginos/neutralinos. Their effect can be easily included in the 
threshold terms, as well as in the running of the f3- and 7-functions, where they appear 
as 0{g i ) g 2 g' 2 , g' A ) terms and are thus numerically unimportant. In Ref. [14[] we have 



shown that, in the case of a heavy supersymmetric spectrum, our analytical expressions 
reproduced the Higgs mass spectrum with an error of less than 2-3 GeV. It can be 
easily checked that light charginos and neutralinos can increase the Higgs mass in < 
2-3 GeV [3[], with respect to heavy ones, which is indeed within the errors of our 
different approximations. Nevertheless for completeness, we include in appendix B the 
leading-log 0(g 4 , g' 2 g 2 , g' 4 ) chargino and neutralino contributions to the CP-even Higgs 
masses and the mixing angle. Throughout this work, we have also implicitly assumed 
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that the gluino masses are of order M s . If the gluinos were, instead, much lighter 
than the characteristic squark masses, the running of the third generation Yukawa 
couplings would be different, inducing a small indirect effect on the two-loop Higgs mass 
computation. The third generation Yukawa coupling running is also modified by the 
presence of light charginos/neutralinos in the spectrum. These two loop contributions 
to the CP-even Higgs masses are also presented in appendix B. 

Our present analysis reproduces, with a high level of accuracy, the values of the 
Higgs masses and mixing angles, for the previously studied case of degenerate left- and 
right-handed squark mass parameters, and for values of the squark left-right mixing 
mass parameters lower than the ones giving the maximal values of the lightest CP-even 
Higgs mass. This comparison holds up to a tiny difference coming from the inclusion in 
this work of the small D-term threshold contributions and vacuum-polarization effects. 
This confirms previous results on the upper bound on the lightest Higgs boson mass 
in the MSSM. 
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Appendix A 

In this appendix we provide the analytical expressions for the Higgs-boson self-energies 
in Eq. ( |3.41| ). We first define the different self-energy contributions as 



AII V (MJ) = AIlW(Mj) + AnW(JWj) + An^)(Mj) + AII^M 2 ) 
where ip = h,H,A. The different contributions in (A.l) read: 
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The different coefficients in (A. 8) and (A. 9) are: 
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where the matrices Z}} and Z^ are those diagonalizing the stop and sbottom squared 
mass matrices, ( |1.3| ) and ( |1.4j ), respectively 
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which arises from the integration of the loop of (scalar) particles, is given by 
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Appendix B 

In this appendix we present the leading-logarithmic 0(g 4 , g 2 g' 2 , g' A ) chargino and neu- 
tralino contributions to the CP-even Higgs masses and mixing angle. We first present 
the modification of the quartic couplings, Eq. fl3.4|), obtained while changing the 
chargino and neutralino masses from Msusy to an overall mass m^. We get, 
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where in the case of charginos lighter that the top quark mass one has to replace m 2 ~ 



by m 2 . From Eq. (B.l) the corresponding contributions to the Higgs mass matrix 
elements can be obtained: 
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Since the chargino and neutralino loops lead to small corrections to the Higgs spec- 
trum, it is possible to obtain approximate expressions for the CP-even Higgs masses, 
while keeping only terms linear in the mass matrix corrections AAifj. The corrections 
to the Higgs mass eigenvalues are given by 

H = v [AAi =F (cos 2a cos 2(3 AX 1 + sin 2a sin 2(3A(X 3 

Analogously, the correction to the Higgs mixing angle a reads 
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In the above, we have always assumed that the corrections to the Higgs mass matrix 
elements are small in comparison not only to the squared Higgs mass eigenvalues but 
also to their difference, m 2 H — m\. If these conditions are not fulfilled, the above 
approximate expressions, Eqs. (B.3) and (B.4) are no longer valid. This might happen, 
for example, in the large tan/3 regime, where for some particular value of the CP- 
odd mass, the CP-even Higgs bosons become almost degenerate. In this case the full 
diagonalization of the Higgs mass matrix is needed in order to get the correct expression 
for the Higgs mass eigenvalues and mixing angle. 

The presence of light gauginos and Higgsinos can also affect the two loop corrections 
to the Higgs masses. This is due to the fact that they affect the top and bottom 
Yukawa coupling beta functions, which enter into the masses through the derivative of 
the quartic coupling beta functions, fy, Eq. ( |3.4|) , and through the soft supersymmetry 
breaking and fermion mass dependence in Eq. ( |3.9|) . The most relevant effects come 
from the possible presence of a gluino lighter than the characteristic stop and sbottom 



masses, m~ < Ms- 



The dominant corrections to (3[ are given by 
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where X, Y = Q, U, D, g, x, t, the superscript XY indicating that we are considering the 
behaviour of the quartic coupling between the energy scales mx and my (mx > my) 
and the value of the functions 9{ should be set to one if m, < mx and 9{ = otherwise. 
The Qt factor is included to make explicit the fact that we are evaluating all quartic 
couplings at the scale m t . The two-loop corrections to the Higgs mass matrix elements 
coming from the variation of the quartic couplings at m t for the case of light gauginos 
and/or Higgsinos may be obtained through Eq. ( |B.2| ), where the variation of the 
quartic couplings may be obtained by a simple generalization of Eq. ( |3.4j ) to the case 
under consideration: 

where XY denote pair of masses in hierarchical order and ti = log(mf/m|). 
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Figure 1: Plot of the Higgs mass mh from Eq. ( |2.2S|) (solid lines), the RG improved 
one-loop MSSM effective potential (dashed lines) and the MSSM effective potential 
considered at Q 2 = m 2 (dotted lines), as described in section 2, for M t = 175 GeV, 
rriA = Ms, and A t = fi = 0. The lower set corresponds to tan/3 = 1.6, and the upper 
set to tan (3 = 15. 
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Figure 2: Plot of the Higgs mass rrih as a function of A t , from Eq. ( |2.29| ), using the 
exact threshold function (|2.21|) (solid line) and the approximation ( p. 22] ) (dashed line), 
for M t = 175 GeV, \i = and rriA = Ms = 1 TeV. The upper set corresponds to 
tan/3 = 15 and the lower set to tan/3 = 1.6. 
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Figure 3: Plot of the Higgs mass m^ as a function of itlq (solid line), for M t = 175 
GeV, tan/5 = 1.6, fi = Ab = 0, my = mo = 1 TeV, and values of the CP-odd Higgs 
mass mA and the stop mixing mass parameter A t equal to mg. The dashed and dotted 
lines denote the same as in Fig. I. 
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Figure 5: Plot of the pole Higgs mass as a function of toq, for M t = 175 GeV, 
tan/3 = 1.6, \i — Aj, — 0, my = itld = 1000 GeV and itla = 300 GeV. The different 
lines denote different values of the A t parameter. Starting from below at tjiq = 1 TeV, 
A t = 0, 0.6, 1.0, 1.5, 1.8, 2.0 and 2.4 TeV, respectively. 
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Figure 6: The same as Fig. 5 but for m\j = 500 GeV and A t = 0, 0.25, 0.5, 0.75, 1.0, 
1.25 and 1.5 TeV, respectively. 
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Figure 7: The same as Fig. 5 but for m v = 100 GeV and A t = 0, 0.2, 0.4, 0.6, 0.8, 
and 1 TeV, respectively. 
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Figure 8: Plot of the pole (solid lines) and running (dotted lines) Higgs masses for a 
running CP-odd Higgs mass itia = 100 GeV, tan (3 = 1.6, tuq = 365 GeV, rrtu = m D = 
500 GeV, Ab = fi = and M t = 175 GeV. The lower (upper) set corresponds to the 
scalar Higgs boson h {H), and middle curve to the pseudoscalar Higgs boson A. The 
mass of the lightest stop is plotted as the dashed line. 
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Figure 9: Plot of the pole Higgs mass M h as a function of mo, for M t = 175 GeV, 
tan/5 = 60, A b = 0, m Q = m v = 1 TeV, A t = 0, 1.5, 2.4 TeV (from bottom to top) 
and /j, — 1 TeV (solid curves) and /x = 2 TeV (dashed curves). 
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